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We study the decays of the lightest supersymmetric particle (LSP) in models with spontaneously 
broken R-parity. We focus on the two cases that the LSP is either a bino or a neutral singlet 
lepton. We work out the most important phenomenological differences between these two scenarios 
and discuss also how they might be distinguished from explicit R-Parity breaking models. In both 
cases we find that certain ratios of decay branching ratios are correlated with either the solar or the 
atmospheric (and reactor) neutrino angle. The hypothesis that spontaneous R-Parity violation is 
the source of the observed neutrino masses is therefore potentially testable at the LHC. 



I. INTRODUCTION 

Experiments on solar and atmospheric neutrinos have demonstrated that neutrinos have non-zero masses 
and mixing angles P, 0, Q . Indeed, with the latest experimental data by the MINOSH and KamLAND Q 
collaborations neutrino oscillation experiments truly have entered the precision era [6j. 

Many models have been proposed, which potentially can explain these data. Certainly the most popular 
are variants based on the see-saw mechanism 0, H S E El- However, the classical seesaw mechanism 
puts the scale of lepton number violation near the grand unification scale and, therefore, can not be directly 
tested. Indirect tests of the seesaw might be possible, if low-scale supersymmetry is realized in nature 
[l2l [Til Hil . H5I [l6| . Measurements of masses and branching ratios of supersymmetric particles can be used 
to obtain indirect information on the r ang e of allowed seesaw parameters, if some specific scenarios for the 
soft-breaking parameters are assumed [l3[ [3, [H[ . On the other hand, a number of neutrino mass models 
exist, in which neutrino masses are generated at the electro- weak scale. Examples include the Babu-Zee 
model [TtJ , Leptoquarks or in case of supersymmetry the breaking of R-parity . These low energy 
models usually lead to testable predictions for future colliders. For example, in the case of R-parity violation 
), in particular for models based on bilinear R-parity breaking terms [20(, correlations between certain 
branching ratios of the lightest sup ersymmetric particle decays and the measured neutrino mixing angles 
have been found [U [H, M, HI • 

In _^pthe neutralino as a dark matter candidate is lost. Recent WMAP data (26|, however, have confirmed 
the existence of non-baryonic dark matter and measured its contribution to the energy budget of the universe 
with unprecedented accuracy. Thus, in Jjt p one needs a non-standard explanation of dark matter (DM). 
Examples for DM candidates in Ifip include (i) light gravitinos [13, [H, [2t| , (ii) the axion [3(| Hl[ or (iii) its 
superpartner, the axino [32l. [33|. to mention a few. 

Whether R-parity is conserved or not can, in principle, be easily decided in the case of explicit R-parity 
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violation since (a) neutrino physics implies that the lightest supersymmetric particle (LSP) will decay inside 
a typical detector of existing and future high energy experiments |2~i1 . [25j and (b) the branching fraction for 
(completely) invisible LSP decays is at most 0(10%) and typically smaller [2ll [34| . Spontaneous violation of 
R-parity (s-lfip ) [35L [36j. on the other hand, implies the existence of a Goldstone boson, the majoron (J). In 
s-$ p the lightest neutralino can then decay according to \° ~^ J +v, i-e. completely invisible. It has been 
shown [34| that this decay mode can in fact be the dominant one, with branching ratios close to 100 %, 
despite the smallness of neutrino masses, in case the scale of R-parity breaking is relatively low. In this limit, 
the accelerator phenomenology of models with spontaneous violation of R-parity can resemble the MSSM 
with conserved R-parity and large statistics might be necessary before it can be established that R-parity 
indeed is broken. 

In this paper we study scenarios where the lightest neutralino is the LSP focusing on two representative 
examples: (i) a bino-like LSP as this case has been extensively studied in the literature, both in case of 
conserved R-parity as well as (explicitly) broken R-parity. We will re-iterate our previous result [34] that 
there are regions in parameter space, where it will be difficult to obtain clarity about the underlying model 
in the first years of LHC. We will also discuss, how measurements of branching ratios can lead to tests of 
the model as the origin of the neutrino mass, in case sufficient statistics for the final states with charged 
leptons can be obtained, (ii) A singlino-like LSP. This case has not been studied in the past but it is the 
only part of the allowed parameter space where singlinos can be produced and studied at future accelerators. 
In addition, this case allows interesting cross-checks with respect to neutrino physics different from a bino 
LSP. 

The paper is organized as follows: in Sect. [II] we present the model with special emphasis on neutrino 
physics as well as the most important couplings of the lightest neutralino. In Sect. InTl we discuss production 
and decays of the lightest neutralino at the LHC as well as a future International Linear Collider (ILC) 
putting some emphasis on how to distinguish bino and singlino LSPs. In Sect. IIVI we discuss the different 
correlations between ratios of neutralino decay branching ratios and neutrino mixing angles. Finally, we 
draw our conclusions in Sect. |Vj 

II. SPONTANEOUS R-PARITY VIOLATION 
A. Model basics 

Spontaneous breaking of a global symmetry leads to a Goldstone boson, in case of lepton number breaking 
usually called the Majoron. Spontaneous breaking of R-parity through left sneutrinos [35j . produces a doublet 
Majoron, which is ruled out by LEP and astrophysical data [UOSl. To construct a phenomenologically 
consistent version of s-lfip it is therefore necessary to extend the particle content of the MSSM by at least 
one singlet field, 9 C , which carries lepton number. For reasons to be explained in more detail below, the 
model we consider (36j j contains three additional singlet superfields, namely, P c , S and $, with lepton number 
assignments of L = —1,1,0 respectively. 

The superpotential can be written as [36| 

W = ^Q i %H u + h%Q i D j H d + h%L i EjH d 

+ h%V c H u -h H d H u $ + h$>V c S+^ 3 . (1) 

The basic guiding principle in the construction of Eq. |T]) is that lepton number is conserved at the level 
of the superpotential. The first three terms are the usual MSSM Yukawa terms. The terms coupling the 
lepton doublets to v c fix lepton number. The coupling of the field <J> with the Higgs doublets generates an 
effective /i-tcrm a la Next to Minimal Supersymmetric Standard Model [38|. The last two terms, involving 
only singlet fields, give mass to D c , S and $, once $ develops a vacuum expectation value (vev). 

For simplicity we consider only one generation of v c and S. Adding more generations of v c or S does 
not add any qualitatively new features to the model. Note also, that the superpotential, Eq. ((T|), does 
not contain any terms with dimension of mass, thus potentially offering a solution to the /i-problem of 
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supersymmetry. The inclusion of S allows to generate a "Dirac"-like mass term for P c , once <& gets a vev. 
The soft supersymmetry breaking terms of this model can be found in [3!| . 

At low energy, i.e. after electro- weak symmetry breaking, various fields acquire vevs. Besides the usual 
MSSM Higgs boson vevs Vd and v u , these are ($) = (D c ) = vr/^/2, (S) — vs/V2 and (z>j) = 

Note, that vr 7^ generates effective bilinear terms = h l v VR/\/2 and that vr, vs and VL t violate lepton 
number as well as R-parity. 



In the basis 



B. Neutralino- neutrino mass matrix 



-i\ 3 , H d ,H u ,is e , v T ,v c , S, <l) 



(2) 



the mass matrix of the neutral fermions following from Eq. (JTJ) can be written as [3£| |4C 
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111 



X" 
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va. x o v 
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M s * 







M s * 


M # 



(•3) 



Eq. ([3]) can be diagonalized in the standard way, 

M N = ATMnAT 1 . 



(4) 



We have chosen the basis in eq. such that N reduces to the MSSM neutralino rotation matrix in the 
limit where (a) R-parity is conserved and (b) the field <J> is decoupled. The various sub-blocks in eq. ([3]) are 
defined as follows. The matrix M % o is the standard MSSM neutralino mass matrix: 



M x o = 





M 2 



/Mi 


- \9' v d +\gvd 
\g'v u -\gv u -n 



\g'v d +±"' 

\gvd - 



g Vu \ 
2 g v u 

A* 

/ 



(5) 



Here, [i = hov^/\^2. m x o ly is the R-parity violating neutrino-neutralino mixing part, which also appears in 
explicit bilinear R-parity breaking models: 



T 

m x <V 



/ -\g'v Le \gv Le e e \ 

-\g'vLu \gvhv e M 
\ -\g'v L T \gvLr e T ) 



(6) 



where vl% are the vevs of the left-sneutrinos. 



in , 



v = fo, 0, 0, -j^Y^Kvu 



(7) 
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and m^o* is 



i£o* = (0 > 0,--j=hov u ,--j=h Vd). (8) 



The "Dirac" mass matrix is defined in the usual way: 



And, finally, 



(m D )< = \K v u- (9) 



(M„c S ) = -l=hv<!>, M u c& = -^=hvs, M s # = -^=hv R , M* = -^=v$. (10) 



The matrix, eq. ([3]), produces ten eigenvalues with vastly different masses. First, since Det(M-yi) = 0, one 
state is massless at tree-level. Then there are two more very light states, together they form to a good 
approximation the three observed, light doublet neutrinos. Their masses and mixing will be discussed in 
detail in the next subsection. 

The remaining seven eigenstates are typically heavy. They can be sub-divided into two groups: Mainly 
doublet and mainly singlet states. There are usually four states which are very similar to the well-known 
MSSM neutralinos. Unless ho is large and Xvg> small, mixing between the phino and the higgsinos is small 
[57[, and there are three singlets. From these singlets, unless (ft- — A) < vr/v$,, v c and S form a quasi-Dirac 
pair, which we will loosely call "the singlino" , Si t 2 — ^75 (^ T ■ Note, that this is a different state compared 

to the NMSSM singlino [4l| which corresponds to <5 in our notation. 

Which of the seven, heavy states is the lightest depends on a number of unknown parameters and can not 
be predicted. In our analysis below we will concentrate on two cases: (a) As in mSugra motivated scenarios 
Mi is the smallest parameter and the lightest state mainly a bino. We study this case in order to work out 
the differences toft) the well-studied phenomenology of the MSSM; and (ii) to the explicit R-parity violating 
case studied in (21J. The second case we consider is (b) the singlino 5 being the lightest state. This case 
is interesting, since it is the only part of the parameter space, where singlets indeed can be produced and 
studied at accelerators. 



C. Neutrino masses 



Since neutrino masses are much smaller than all other fermion mass terms, one can find the effective neu- 
trino mass matrix in a seesaw-type approximation (39l . |40| . First we define the small expansion parameters 
£y, which characterize the mixing between the neutrino sector and the seven heavy neutral fermion states, 
the "neutralinos" of the model, 



£ = m 3X 7-M H . 

The sub-matrix describing the seven heavy states of eq. © is 



M„c S 
M„c S M s * 



(11) 



/ M x o 



M H = 



M^c # M s * 



Ma 



(12) 



and 



m 3X 7 = ( m £v m D 



(13) 
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We have neglected m x o„c in eq. (fT2|) since it is doubly suppressed. The "effective" (3,3) neutrino mass 
matrix is then given in seesaw approximation by 



(70 = -m 3X 7-MH m^ x7 . 



(14) 



In the following we will use the symbol Nij with i,j = 1, ...7 as the matrix which diagonalizes eq. (fT2| . Our 
N reduces to the MSSM neutralino mixing matrix N, in the limit where the singlets decouple, i.e. h —> 
or M$ — > oo. After some straightforward algebra £y can be written as 



where the effective bilinear R-parity violating parameters and A, are 

1 K V2 



and 



The coefficients K are given as 



x 2g'M 2f x 

K A = ™ ° 



Ki = 



2gMifi 



-v u a + 



2v u 



K A = v d a+- 

K 5 - VRb 
~2vu 

6 v sb 



K 



A 



2v u 
v<s,b 
2^ 



K\ 
K? 

k! 

Kf 
K\ 

Kt 
Kl 



A t = tiV d + [iv Lt . 
2g'M 2f M 1 



(15) 
(16) 

(17) 

(18) 



Illy 

2gMm 



c ADet(M H )a s 
( 7^ 



h 2 [iv 



2v u 



4Det(M H 
vrc 
2v u 

c 

2^/2v u v R vq>h h 2 
v$c 
2v u 



■(AM 1 M 2 fiv u - m 7 v d v 2 ) 



[ + V2hv^v R v s 



2y/2Det(M H )b 2 



ft 3 



m 7 v u v R v<f> 



The coefficients a, b and c are defined as 



4V2Det(M H ) { - hVRVS+ 2 Xvl + h ° VdVu) > 



(19) 



b = 



m 1 h 2 \i 
4Det(M H ) 



Vu{v 2 u -v 2 d ), 



h 2 n 



-v 2 l (2M 1 M 2 fJ, - m 1 v d v u ). 



Det(M H ) 

Det(Mjj) is the determinant of the (7, 7) matrix of the heavy neutral states, 

Det(M H ) = -^h h 2 vl [A(2M 1 M 2 fi - m 1 v d v u )(-hv R v s + ^Au| + h v d v u ) - /i m 7 (u 2 - v 2 d ) 2 ] 



(20) 



and v 2 = v 2 a + v d . The "photino" mass parameter is defined as m 7 = g 2 M\ + (f 2 M 2 . Note that the K\ and 
K\ reduce to the expressions of the explicit bilinear R-parity breaking model | 42] , in the limit M$ — ► oo and 
in the limit h, ho — > 0, i.e. b = c = 0. 
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The effective neutrino mass matrix at tree-ievel can then be cast into a very simple form 

- i me Jl)i] = aA * A j + + £ 3 K i) + ce i e i ■ ( 21 ) 

Equation (|21[) resembles very closely the corresponding expression for the explicit bilinear R-pa rity breaking 
model, once the tree-level and the dominant 1-loop contributions are taken into account [12, [H, EH ■ Eq. 
pT|) reduces to the tree-level expression of the explicit model [58| 

< m "»« - iijSv AiAi (22) 

in the limit Af$ — > oo and in the limit h, ho — ► 0. Different from the explicit model, however, the spontaneous 
model has in general two non-zero neutrino masses at tree-level. With the lightest neutrino mass zero at 
tree-level, the s-^ p model could generate degenerate neutrinos only in regions of parameter space where 
the two tree- level neutrino masses of eq. (|2Tj) are highly fine-tuned against the loop corrections. We will 
disregard this possibility in the following. 

Neutrino physics puts a number of constraints on the parameters Aj and e$. However, in the spontaneous 
model there is no a priori reason which of the terms gives the dominant contribution to the neutrino mass 
matrix, thus two possibilities to fit the neutrino data exist: 

• case (cl) A generates the atmospheric mass scale, e*the solar mass scale 

• case (c2) e generates the atmospheric mass scale, A the solar mass scale 

The absolute scale of neutrino mass requires both |A|//n and |e|//i to be small, the exact numbers depending 
on many unknown parameters. For typical SUSY masses order 0(100 GeV), |A|//i 2 ~ 10~ 6 - 10~ 5 . If some 
of the singlet fields are light, i.e. have masses in the range of 0(0.1 — few) TeV, also |ei//ii| can be as small 
as \e\/ fJ, ~ 10~ 6 -10~ 5 . On the other extreme, independent of the singlet spectrum, |e|//x can not be larger 
than, say, |e|//i ~ 10~ 3 , due to contributions from sbottom and stau loops to the neutrino mass matrix 

The observed mixing angles in the neutrino sector then require certain ratios for the parameters Aj/Aj and 
Ci/ej. This can be most easily understood as follows. As first observed in (45[, the so-called tri-bimaximal 
mixing pattern 



t/ HPS 



/ ,fl ,fl o 

V 3 V 3 u 



V 



i i i_ 

\/6 \/3 V2 



(23) 



V6 4 V2 

is a good first-order approximation to the observed neutrino angles. In case of hierarchical neutrinos M.^ a9 — 
(0,m, M), where m (M) stands for the solar (atmospheric) mass scale, rotating with C/ HPS to the flavour 
basis leads to the following neutrino mass matrix 

i/0 \ -t / m m m\ 
m hps = - M —M + - to m m . (24) 
2 y —M M ) 3 y to m to J 

In case the coefficient b in eq. (|19|) is exactly zero, i.e. for tan /3 — 1, the model would produce a tri-bimaximal 
mixing pattern for Ai = 0, A2 = —A3 and t\ — (.2 = £3, in case (i). For case (ii) the conditions on A^ should 
be exchanged with the conditions for the &j and vice versa. 

In reality, since tan/3 ^ 1 in general, neither is b exactly zero, nor need the neutrino mixing angles be 
exactly those of eq. (|2"3"]) . One then finds certain allowed ranges for ratios of the A,; and ej. In case (i) one 
gets approximately 

-tan 2 6» R , (25) 



Aj 



VAl+A 3 



2 



" x ~'tan 2 6» A tm, 



tan 2 1 
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Here, e = Uj ■ e*with {U U ) T being the matrix which diagonalizes the (3,3) effective neutrino mass matrix. 
In case (i) (U V ) T is (very) approximately given by 



|A| 





A1A2 



Ai A 3 



V 



|A| 



VAf+A||A| 
A 3 

VAf+Af 
A2 

|A| 



A2 

VaI+a| 
A3 

|A| 



(26) 



Note that U][ is the matrix which diagonalizes only the part of the effective neutrino mass matrix proportional 
to AjAj. Again, for the case (ii) replace A, <-> in all expressions. 



D. Approximated couplings 

With R-parity violated the lightest supersymmetric particle decays. Here we list the most important 
couplings of the lightest neutralino in the seesaw approximation. In the numerical calculation discussed 
in the next sections, we always diagonalize all mass matrices exactly and obtain the exact couplings. For 
the understanding of the main qualitative features of the LSP decays, however, the approximated couplings 
listed below will be very helpful. 

We define the "rotated" quantities: 

Xi = (U v ) T ik x k , yn = (U u ) T ik y k:j . (27) 
X? — W ± — If couplings are found from the general expressions for the x° — W ± — x T vertices 

£ = Xir (0%?P L + 0%?P R ) X °jW- + xh» (O n dfP L + 0%gP a ) X J W+ (28) 

as 

~ cnw 1 , v d rgv d N 12 + M 2 N 14 g{2^ 2 +g 2 v u v d )N 12 +g 2 v u (fi + M 2 )N u , 
0«i = 2 9{kEh ^ Jb 61 + 2^7 A ^ 

Dei + is the determinant of the MSSM chargino mass matrix. Here, 

ffil — ^ffil J ■ 

The Lagrangian for Xi ~ X.'j ~ Z 

c = Ixh^oi^PL + O n R l z P R ) X ° 3 Z, (31) 

gives for Xi — v % — Z 

01% = - — a [iaNu + i a N 12 + 2£ H N U + £ i5 N 15 + £ i6 N w + £„N ir ] , (32) 
2 cos t>w 

Onnz ( r~}nnz\* 

Ril — -\ u Lil) ■ 

The most important difference to the explicit R-parity violating models comes from the coupling Xi~X*j — P k 

£ = \tipZ%PL + O^Pn^Pl (33) 
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with 

OIZj - R'jmO'ZL (34) 

U Ri\J — \ U Lilj) ■ 

Because the spontaneous breaking of lepton number produces a massless pseudo-scalar, eq. (j3"3")l leads to a 
coupling Xi ~ J ~ v i, i-e a new invisible decay channel for the lightest neutralino. For vl <C vr, v$ one can 
find an approximation to the Majoron which, in leading order, is given by 

(0,0,^,0,^,-^). (35) 



Here, V = \Jv\ + w| and terms of order have been neglected. The "unrotated" couplings O^"™^ are 

2 



°'lTz° = - l M ki (9'Nn- 9 N 12 ), (36) 



°l*s = ^h(i i5 N 17 + i i7 N m ), 

OlZ* = -i—Nu + ^=h(£ i6 N 17 +i i7 N m ). 
V R v 2 

In the limit vr, vs — > oo one can derive a very simple approximation formula for O^o uk j. It s given by [Bl^ ] 

\O nVhJ \ ~ -^iV 14 + ^(ff'iVu - ff JV u ) + ft.O. (37) 

Eq. (|3T[) serves to show that for constant e and «l, O^o„ fc j — > as vr goes to infinity. This is as expected, 
since for vr — ► oo the spontaneous model approaches the explicit bilinear model. Note, that only the presence 
of the field v c is essential for the coupling Eq. IpT?]) . If 5 is absent, replace V — > ur. 

In addition to the Majoron in considerable parts of the parameter space one also finds a rather light singlet 
scalar, called the "scalar partner" of the Majoron in [4(J, Sj. From the Lagrangian 

£ = \t{oz] k p L + ()■;■;>, I'uj^s?. (38) 

one finds the coupling x\ — Sj ~ v i as 

Onns r>s n'nns /"QfA 

LilSj ~ n Sjk u Lilki 

rynns / r\nns \ * 

^miSj — v J LiiSj) ■ 

Different from the Majoron, however, there is no simple analytical approximation for Rsj- We write sym- 
bolically 

RS Sjk = {Rs J H d ,RsjH u ,R Sj Lo,Rsj<S>,R S jS' R S.ji>' : )> (40) 

and define unrotated couplings by 

O'^Sh = \ [(9^2 - 9'ia)N 13 + (gN 12 - g'N n )£ i3 - V2h (£ i7 N u + £ i4 N 17 )] , (41) 

O^ik = J [(sfiii - git2)N u + (g'Nu - gN 12 )U - V2h Q (£ z7 N 13 + £ i3 N 17 ) - V2^N 15 ] , 
2 vr 

Ol"i4 = ^= [ - hoGisNu, + iuNvt) + h(£ m N 15 + £ i5 N 16 ) + \£ i7 N 17 ] , 

°l*s = ^h(£ i5 N 17 + t„N 15 ), 

<fe = -—N u + ^h(£ i6 N 17 + £ i7 N 16 ). 
vr V2 
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As eqs fill) shows, \i ~ * $J + v i has a partial decay width similar in size to the decay xi - * + v i, as soon 
as kinematically allowed. Since, on the other hand, Sj decays practically always with a branching ratio close 
to 100 % into two Majorons, Xi ~* Sj + v% gives in general a sizeable contribution to the invisible width of 
the neutralino. 

Finally, we give also the coupling X° — J — Xj> f° r the case of two heavy neutralinos. Here, 

Ouh = -^[MN j7 N i5 + N„N j5 ) ~ v R (N j7 N i6 + N, l7 N j6 )} . (42) 



III. LSP PRODUCTION AND DECAYS 

In this section we discuss the phenomenology of a neutralino LSP in s-I/Lp at future colliders. We do not 
attempt to do an exhaustive study of the (quite large) parameter space of the model. Instead we will focus 
on the most important qualitative differences between s-Ijip , the previously studied case of explicit bilinear 
Ijtp [2l|, [4^, [43, S3] and the MSSM. All numerical results shown below have been obtained using the program 
package SPheno [4!|, extended to include the new singlet superfields v c , S and $. 

Unless mentioned otherwise, we have always chosen the JjL v parameters in such a way that solar and 
atmospheric neutrino data Q are fitted in the correct way. The numerical procedure to fit neutrino masses 
is the following. Compared to the MSSM we have a number of new parameters. For the superpotential of eq. 
([T]) these are ho, h and A, as well as the neutrino Yukawas h l v . In addition, there are in principle also the soft 
SUSY breaking terms, which generate non-zero vevs, Vr, vs, v$ and vj Ji for v c , S, $ and i>i, respectively. 
We trade the unknown soft parameters for the vevs. For any random choice of MSSM parameters, we can 
reproduce the "correct" MSSM value of /i for a random value of i>$, by appropriate choice of ho- For any 
random set of h, A, vs and vr, we can then calculate those values of h\, and u^, using eq. (j2Tj) . such that 
the corresponding ti and Aj give correct neutrino masses and mixing angles. There are two options, how 
neutrino data can be fitted, i.e. the cases (cl) and (c2), defined in section Hi CI We discuss the differences 
between these two possibilities below. 

In the following we will study only two 'limiting' cases, which we consider to be the simplest possibilities 
to realize within the parameter space of the model: (a) a bino-like LSP and (b) a singlino LSP. We note, 
however, that theoretically also other possibilities exist at least in some limited parts of parameter space. 
For example, one could also have that the phino, <&, is the lightest R p odd state. However, with the 
superpotential of eq. (JTJ) , for any given value of \i, v$ has a minimum value. Since the product Av$ also 
determines approximately the phino mass, a very light phino requires a certain hierarchy A <C ho, h, which 
might be considered to be a rather special case. Also in mSugra in the region where mo is large one can find 
points in which /1 ~ Mi and the lightest (MSSM) neutralino has a significant higgsino component. Since 
both, a higgsino as well as a phino LSP show some differences in phenomenology compared to the bino and 
singlino LSPs discussed here, we plan to study higgsino and phino LSPs in a future publication. 

A. Production 

Since neutrino physics requires that the R-parity violating parameters are small, supersymmetric produc- 
tion cross sections are very similar to the corresponding MSSM values, see for example [50| and references 
therein. Over most of the MSSM parameter space one expects that mainly gluinos and squarks are directly 
produced at the LHC and that the lightest neutralinos appear as the "final" decay products at the end of 
possibly long decay chains of sparticles. In addition charginos, neutralinos and sleptons can be produced 
directly via Drell-Yan processes provided that they are relatively light. 

Cross sections for direct production of singlinos are always negligible. There are essentially two possibilities 
how singlinos can be produced in cascade decays. Firstly, a somewhat exotic chance to produce singlinos 
occurs if at least one of the MSSM Higgsinos is heavier than <!> and both ho and h are large. In this case 
Si appear in decay chains such as H u ^d — > $ + A"i — > S + X2, where Xi denotes the additionally produced 
particles. Secondly, there is the possiblity that singlinos are the LSPs. Squarks and gluinos will then decay 
fast to the NLSP, which then decays to S. A typical decay chain might beq^q + B^q + 6>i,2 + J. Other 
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NLSPs such as, f\ will decay mainly via t\ — > S1.2 + t, i.e. again ending up in singlinos. The total number 
of singlino events therefore will be simply approximately equal to the number of SUSY events for singlino 
LSPs. 



B. Decays 




90 100 110 120 90 100 110 120 

hv<s,/V2 [GeV] [GeV] 



Figure 1: Masses of the three lightest neutralinos (left) and branching ratios for the most important decay modes of 
the lightest state (right) versus -^=hvi£ for a specific, but typical example point. The MSSM parameters have been 
adjusted such that the sparticle spectrum of the standard point SPSla' is approximately reproduced. The singlet 
parameters have been chosen randomly, vr — vs = 1 TeV, e and A have been fitted to neutrino data, such that A 
generates the atmospheric scale and e the solar scale. For a detailed discussion see text. 

Here we will discuss the main decay modes of bino and singlino LSPs. We will first discuss the parameter 
range, where m x o > m w ± , such that two-body decays of \i to gauge bosons are kinematically allowed. Fig. 
(JTJ) shows an example of the three lightest neutralino mass eigenvalues (left) and the main decay modes of 
Xi (right) as a function for fixed values of all other parameters. This point has been constructed in such 
a way, that the MSSM part of the spectrum, all production cross sections and all decay branching ratios, 
apart from the lightest neutralino decays, match very closely the mSugra standard point SPSla' [50] - Here, 
vr — vg = 1 TeV has been chosen as an arbitrary, but typical example. 

The left part of fig. (TTJ shows how the quasi-Dirac pair 6>i,2 evolves as a function of ^i. For low values 

(i.e. < Mi)of this parameter combination Si is the LSP, for large values a B is the LSP. The right side of 
the figure shows the final states with the largest branching ratios. For low values of the LSP mass, J/Sj + v 
is usually the most important, i.e. there is a sizeable decay to invisible final states, even for a relatively high 
vr, see also the discussion for fig. Next in importance are the final states involving W ± and charged 
leptons. Note, that the model predicts 

y\j?r(x;^° + ^) „ 9 U3) 
2 Ei Br (Xi -> W+ + ll) 4cos 2 6V 1 ' 

with g being a phase space correction factor, with g — > 1 in the limit m x o — > 00 [29j . Equation (143|) can 
be understood with the help of the approximative couplings eq. (|29p and eq. (|32|) . The relative size of the 
branching ratios for the final states W + e, W + fj, and W + t depends on both, (a) the nature of the LSP 
and (b) the fit to the neutrino data. We will discuss this important feature in more detail in section ITVl 
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Generally, for m x o > m w ± three-body final states of the neutralino decay are less important than the 
two-body decays shown in fig. (CD). Especially one expects that the final state vbb has a smaller branching 
than in the case of explicit Jjt p [21(. This is essentially due to the fact, that \e\/n is smaller in s-$m with 
a "light" singlet spectrum than in a model with explicit bilinear IjLp , see the discussion in section III CI A 
smaller |e|//u leads to smaller couplings between y^ — l — I, Xi — <7 — Q and especially Xi ~ v ~ h°, see also 
couplings in [2l|. We have checked numerically, that Br(xi — > u + h°), if kinematically open, is typically 
below 1% for singlets in the O(TeV) range. 
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Figure 2: Branching ratios for the most important decay modes of the lightest neutralino state versus -^hv<s> for a 
specific, but typical example point. The MSSM parameters have been adjusted such that the sparticle spectrum of the 
standard point SU4 is approximately reproduced. The singlet parameters have been chosen randomly, vr — vs — 1 
TeV, e and A have been fitted to neutrino data, such that A generates the atmospheric scale and e the solar scale. 
The different final states are as follows. In the left figure, as ordered on the right side, from top to bottom the lines 
are Br(x? — * [invisible]) (full line, blue), Br(x? — > A 1 *??') (short-dashed, red), Br(x? — ► Tqq') (large-dashed, light blue), 
Br(x? — > vqq) (full, yellow), Br(x? — * vbb) (full, pink) and Br(x? — * eqq') (full, green). In the right figure, purely 
leptonic modes, from top to bottom (on the right side): Br(x? — > vfir) (full, yellow), Br(x? — * ve/i) (dashed, green), 
Br(x? -» ver) (full, red), Br(x? -> v\i\i) (dashed, light blue), Br(x? — > vtt) (full, pink) and Br(x? — > ^ee) (dashed, 
darker blue). For a detailed discussion see text. 

For the case of m x o < m^± fig. ([2]) shows an example for the most important final states of the lightest 
neutralino decay as a function of ^=-. As in the fig. (CL| to the left the lightest neutralino is a singlino, to the 
right of the "transition" region the lightest neutralino is a bino. Note that the point SU4 [HTj] produces a bino 
mass of approximately rrig ~ 60 GeV, thus the only two body decay modes which are kinematically allowed 
are J + v and - very often, but not always - Sj + v. One observes that these invisible decay modes have 
typically a larger branching ratio than in the case m x o > m w ± shown in fig. (CD) . This fact is essentially due 
to the propagator and phase space suppression factors for three body decays. For a bino LSP the invisible 
decay has the largest branching fraction. Semileptonic modes are next important with typically liqq' being 
larger than vqq. It is interesting to note, that in the purely leptonic decays, lepton flavour violating final 
states such as fir have branching ratios typically as large or larger than the corresponding charged lepton 
flavour diagonal decays (/xp and rr). These large flavour off-diagonal decays can be traced to the fact 
that neutrino physics requires two large mixing angles. The branching ratios shown in fig. |J5J) should be 
understood only as representative examples - not as firm predictions. Especially for the case of a bino LSP, 
the partial width to the final state J + i.e. invisible final state, can vary by several orders of magnitude, 
see the discussion below. The predictions for relative ratios of the different (partially or completely) visible 
final states is much tighter fixed, because these final states correlate with neutrino physics, as we discuss in 
section IIVI 
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If the S is the LSP, a bino NLSP decays dominantly to the singlino plus missing energy, as is shown in fig. 
©. The final state can be either Si + J or Si + 2 J, the latter due to the chain B — ► <S 2 + J —> Si + 2 J, where 
the 2nd step has always a branching fraction very close to 100%. However, a special opportunity arises if 
h is low. In this case ^ Br{x\ — B — > W ± + If) can easily reach several percent and it becomes possible 
to test the model with the bino decays and the singlino decays at the same time. This would allow a much 
more detailed study of the model parameters than for the more "standard" case where only either singlino 
or bino decay visibly. We note that for any fixed value of h, Y^i Br{x% — B ~> W ± + If) depends mostly on 
vr (and to some extend on «$). Low values if vr lead to low ^\ Br(x® — B — > W ± +lf) as we will discuss 
next. 




Figure 3: Sum over Br(x3 — B —* W ± + If) versus h for MSSM parameters resembling the standard point 
SPSla', random values of the singlet parameters and with the condition of Si being the LSP. The dominant decay 
mode for the B in all points is B — * <Si+}2 , with the missing energy due to either J or 2 J emission. For low values 
of h one can have visible decays of the B reaching (20 — 30)%, for h larger than, for say, h — 0.05 B decays to Si 
plus missing energy with nearly 100%. 

Fig. (|4]) shows the sum over all at least partially visible decay modes of the lightest neutralino versus 
vr in GeV, for a set of u$ values u$ = 10-40 TeV for the mSUGRA parameter point (mo = 280 GeV, 
mi/2 — 250 GeV, tan/3 = 10, A = —500 GeV and sgn(/i) = +). This point was constructed to produce 
formally a f2 x o ft, 2 ~ 1 in case of conserved R-parity, much larger than the observed relic DM density [26| . 

The left plot shows the case x? — B, the right plot Xi — F° r B, Br(B — > J + v) very close to 100 % 
are found for low values of vr. This feature is independent of the mSugra parameters, see the correspoding 
figure in (34|. In this case large statistics becomes necessary to find the rare visible neutralino decays, which 
prove that R-parity is broken. The inconsistency between the calculated £l x oh 2 and the measured SlcDMh 
might give a first indication for a non-standard SUSY model. 

Figure (d]) to the right shows that the case Xi — S has a very different dependence on vr. We have 
checked that this feature is independent of the mSugra point. For other choices of mSugra parameters larger 
branching ratios for Br(5 — > J + v) can be obtained, but contrary to the bino LSP case, the sum over the 
invisible decay branching ratios never approaches 100 %. 

Figure ([5]) shows the calculated decay lengths for the lightest neutralino for the same choice of parameters 
as shown in fig. (|4|). To the left the case \i — B, to the right Xi — <5. Decay lengths depend strongly on 
Vr. Singlinos tend to have larger decay lengths than binos for the same choice of parameters. However, a 
measurement of the decay length alone is not sufficient to decide whether the LSP is a singlino or a bino. If 
the nature of the LSP is known, observing a finite decay length allows a rough estimate of the scale vr, or 
at least to establish a rough lower limit on dr. 

Summarizing this discussion, it can be claimed that observing a decay branching ratio of the LSP into 
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Figure 4: Sum over all at least partially visible decay modes of the lightest neutralino versus vr in GeV, for a set 
of w$ values = 10-40 TeV for the mSUGRA parameter point m = 280 GeV, m 1/2 = 250 GeV, tan/3 = 10, 
Ao = —500 GeV and sgn(fi) — +. To the left Xi — B\ to the right Xi — S- The plot demonstrates that the branching 
ratio into B — > J + v does depend strongly on the value of vr and to a minor extend on v$. Lowering vr one can get 
branching ratios for the invisible decay of the B very close to 100 %, thus a very MSSM-like phenomenology. The 
right plot demonstrates that such a possibility does not exist in the case of an S LSP. 
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Figure 5: Decay length of the lightest neutralino in meter versus vr. To the left: Bino LSP; to the right: Singlino 
LSP. All parameters have been chosen as in fig. @. 



completely invisible final states larger than Br(x" — ► ^ [invisible] ) > 0.1 is an indication for s-$p . Finding 
Br(x? ~~ * X) [invisible] ) ~ 100% shows furthermore that the \i must be a bino and measuring Br(xi — * 
[visible]) for a bino LSP gives an order-of-magnitude estimate of vr. 
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C. Possible observables to distinguish between Singlino LSP and bino LSP 

Since bino and singlino LSP decays have, in principle, the same final states, simply observing some visible 
decay products of the LSP does not allow to decide the nature of the LSP. In this subsection we will 
schematically discuss some possible measurements, which would allow to check for the LSP nature. 

As shown above, if the singlino is the LSP and the bino the NLSP, one can have that for the bino decays 
to standard model particles compete with the decay to the singlino LSP. If both particles, the bino and the 
singlino LSP have visible decay modes, it is guaranteed that the singlino is the LSP. If the bino decays only 
invisibly to the singlino, a different strategy is called for. We discuss two examples in the following. 

In the following discussion we will replace the neutralino mass eigenstates by the particles which correspond 
to their main content to avoid confusion with indices. At the LHC one will mainly produce squarks and 
gluinos which will decay in general in cascades. A typical example is qh — > qW and the wino decays further 
to a bino LSP as for example: 

W — > e + eT — > e~e + B — > e~ e + [iqq (44) 
W -> e+g- -► e-e + B -> e~e + Jv (45) 

In this case one can measure in principle the neutralino mass from the first decay chain. In the invariant 
momentum spectrum of the e + e~ pair the edge must correspond to this mass. In the case of the singlino 

W -► e + e~ -> e'e + B -> e^e+JS -► e~e + Jnqq (46) 
W -> e+g- -> e-e+B -> e~e + JS -> e~e + JJv (47) 

In this case one has on average more missing energy than for a bino LSP. However, in both cases one can 
study spectra combining the jet stemmin g fr om the squark and the e + e~ pair and obtain information on the 
masses from the so-called edge variables 52]. In addition one can use additional variables like, for example, 
rriT2 [53l HU, [HI to obtain information on the LSP mass. Note, that this variable works also if there are 
additional massless particles involved, although at the expense of available statistics [56j . In addition one 
can obtain the invariant mass of the LSP from the final state fxqq. In the case where the LSP has a decay 
length measurable at the LHC, one can separate the latter decay products from the other particles in the 
event and, thus, reduce considerably the combinatorial problems associated with the correct assignment of 
the jets. In the case of a bino LSP one would find that all the three different measurements yield the same 
mass for the LSP. In the case of a singlino LSP, on the other hand, one would obtain that the LSP mass 
reconstructed from the edge variables does not coincide with the mass reconstructed from the fiqq spectrum. 
This would indicate that there are two different particles involved. (Such a difference might also be visible 
in the txiti variable.) However, in all cases detailed Monte Carlo studies will be necessary to work out the 
required statistics, etc. 

Distinguishing bino and singlino LSPs will become considerably easier at a future international linear 
collider. In e + e _ one can directly produce a bino LSP but not a singlino LSP and, thus, the identification 
of the correct scenario should be fairly straightforward. 

IV. CORRELATIONS BETWEEN LSP DECAYS AND NEUTRINO MIXING ANGLES 

Correlations between LSP decays and neutrino mixing angles depend on the nature of the LSP. Above 
we have discussed some possible measurements which, at least in principle, allow to distinguish bino from 
singlino LSPs. In this subsection we assume that the nature of the LSP is known. 

1. Bino LSP 

We note that the following discussion is valid also if the bino is the NLSP which, as discussed above, 
decays with some final, but probably small percentage to visible final states. 
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In explicit bilinear R-parity violation the coupling of the bino component of the neutralino to gauge bosons 
and leptons is completely dominated by terms proportional to Aj, as has been shown in [21]. Although the 
coefficients for the spontaneous model are more complicated, see the discussion in section III D| generation 
dependence for the coefficients for the coupling Xi ~ W — U appear only in the terms Aj and ej, i.e. K\ and 
K\ are independent of the lepton generation. Numerically one finds than that the terms proportional to Aj 
dominate the \\ ~ W — U coupling for a bino LSP always. This is demonstrated in figs. ([6]) and ([7|). Here we 
have numerically scanned the mSugra parameter space, with random singlet parameters and the additional 
condition that the LSP is a bino. For the left (right) figures we have numerically applied the cut Aft > 0.5 
(Aft > 0.9). 

Fig. © [ 0] shows the ratio of branching ratios Br(B -> W+e) /Bt(B -> W+n) [Br(B -> W+n) /Br(B ->■ 
W + t)] versus (A e /A M ) 2 [(A M /A r ) 2 ]. To establish a correlation between ratios of Aj and the bino decay 
branching ratios, a bino purity of Aft > 0.5 is usually sufficient. The figures demonstrate that the correlations 
get sharper with increasing bino purity. 

We have checked that for neutralinos with mass lower than mw one can use ratios of the decays B — > Uqq' 
for the different U in the same way to perform a measurement of Aj ratios. Plots for this parameter region 
are rather similar to the ones shown for the case B — > kW, although with a somewhat larger dispersion, and 
we therefore do not repeat them here. 

With the measurement of ratios of branching ratios different consistency checks of the model can be 
performed. In case (cl), i.e. A explaining the atmospheric scale, the atmospheric and the reactor angle 
are related to W + I final states, as shown in fig. ((Sj> . Here we show the ratios TZ^ — g^|*uZ^jy-) versus 
tan 2 0Atm (left) and lZ e — Br(x t ^eW) versus sin 2 9b (right) for a bino LSP, for an assumed 

bino-purity of Aft > 0.8. The vertical lines are the 3ct cl. allowed experimental ranges (upper bound), 
horizontal lines the resulting predictions for the two different observables 1Z. Given the current experimental 
data, one expects p^i^'wl in the range [0.4,2.1] and , = Br ^i~* eW ) < 0.06. 

' ^ Br(xi^rW) y/ Br( X "^ l-i-W) 2 +Br(x1^rW) 2 ~ 



Different from fig. ©, in case of (c2), i.e. A explaining the solar scale, the ratio 



BrjxI^eW) 



^Brixl^/iWy+Brixl^rW) 2 

correlates with tan 2 # , as shown in fig. Here, from the 3a cl. allowed range of the solar angle one 

expects to find Brjx^eW) ^ [0.25, 0.851. Finding this ratio experimentally to be larger than 

yf Br(x l ^^W) 2 +Br(x° l ^TW) 2 

the one indicated by the solar data, i.e. , „ Br(xi~*eW) >> 1, rules out the model as the origin 

J ' ^/Br( X ( I^^W) 2 +Br(x'I^TW) 2 ' 6 

of the observed neutrino oscillation data. Similarly a low (high) experimental value for this ratio indicates 
(for a bino LSP) that case (cl) [(c2)] is the correct explanation for the two observed neutrino mass scales. 



2. Singhno LSP 

Different from the bino LSP case, for singlinos coupling to a lepton Zj-W pair terms proportional to 
dominate by far. This is demonstrated in fig. (TIT))) , where we show the ratios gr[^uZ^vt?) (left) versus 

(e e /e^) 2 and g^*uZ^j^j (right) versus (e M /e T ) 2 for a singlino LSP. Note that mixing between singlinos and 
the doublet neutralinos of the model is always very small, unless the singlino is highly degenerate with the 
bino. Consequently singlinos are usually very "pure" singlinos and the correlations of the Zj-W with the 
ratios is very sharp. 

Depending on which case, (cl) or (c2), is chosen to fit the neutrino data, the corresponding ratios of 
branching ratios are then either sensitive to the atmospheric and reactor or the solar angle. This is demon- 
strated in figs. (TTT1) and (TT21) . Here, fig. (fTTI) shows the correlation of lZ e — Br{ Xl -^eW) with 

6 l| V H V ! 6 1 V y/ Brixl^^WY+Brixl^rW) 2 

tan 2 8q for the fit (cl). This result is very similar to the one obtained for the fit (c2) and a bino LSP. For 
this reason the nature of the LSP needs to be known, before one can decide, whether the measurement of a 
ratio of branching ratio is testing (cl) or (c2). 

Figure (fl"2"|) shows the dependence of 7£ M = g^ X oZ^j^ versus tan 2 QAtm (left) and lZ e — 
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Figure 6: Ratio g^zf^ versus (A e /A M ) 2 for a bino LSP. To the left: "Bino-purity" Aft > 0.5, to the right: 
Nil > 0.9. All points with tulsp > mw- 
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Figure 7: Ratio ^l[^hZrw] versus (A^/A T ) 2 for a bino LSP. To the left: "Bino-purity" Nfi > 0.5, to the right: 
Nix > 0-9. 



= versus sin 2 8 R (right) for a singlino LSP, using the neutrino fit (c2). Again 

■J Br{x}— »fiH / ) 2 +Br(xi- >Tff) ! 

one observes that this result is very similar to the one obtained for a bino LSP and fit (cl). This simply 
reflects that fact, that neutrino angles can be either fitted with ratios of or with ratios of Aj and singlinos 
couple mostly proportional to e*, while binos are sensitive to A^. 

In case the singlino is the LSP and the bino, as the NLSP, decays with some measurable branching ratios 
to W — k, both Aj and e$ ratios could be reconstructed, which would allow for a much more comprehensive 
test of the model. 
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for a bino LSP. "Bino-purity" > 0.8. Vertical lines are the 3<r c.l. allowed experimental ranges, horizontal lines 
the resulting predictions for the fit (cl), see text. 




Figure 9: Ratio TZ e = —, Br( Xl -^eW) versus tan 2 6>© for a bino LSP. "Bino-purity" Nu > 0.8. Vertical 

lines are the 3cr c.l. allowed experimental ranges, horizontal lines the resulting predictions for the fit (c2), see text. 



V. CONCLUSIONS 



We have studied the phenomenology of a neutralino LSP in a supersymmetric model in which neutrino 
oscillation data is explained by spontaneous R-parity violation. We have concentrated the discussion on the 
case that the LSP is either a bino, like in a typical mSugra point, or a singlino state, novel to the current 
model. We have worked out the most important phcnomcnological signals of the model and how it might 
be distinguished from the well-studied case of the MSSM, as well as from a model in which the violation of 
R-parity is explicit. 

There are regions in parameter space, where x° decays invisibly with branching ratios close to 100 %, 
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allowed experimental ranges, horizontal lines the resulting predictions for the fit (cl), see text. 



despite the smallness of neutrino masses. In this limit, spontaneous violation of R-parity can resemble the 
MSSM with conserved R-parity at the LHC and experimentalist would have to search for the very rare 
visible decay channels to establish the R-parity indeed is broken. 

The perhaps most important test of the model as the origin of the observed neutrino masses comes from 
measurements of ratios of branching ratios to JU-boson and charged lepton final states. Ratios of these 
decays are always related to measured neutrino angles. If SUSY has a spectrum light enough to be produced 
at the LHC, the spontaneous model of R-parity violation is therefore potentially testable. 



19 







M 


Jmm • 




••JSC 

. •!>•: 

r • . • 

• 







10" 



10 u 



tan 2 ^Atm 



10" 



I 



-2 



Icq 
> 



10 



-3 



•./TO?"' 

• * •...•••*".*/ 
• • .** i • 



10 



-3 



10 



-2 



10" 



sin 2 6^ 



Figure 12: Ratio TZn = i^rir^SQ versus tan 2 0Atm (left) and 7J. e = . „ Br(i "" >eff ) versus sin 2 fl 

(right) for a singlino LSP. Vertical lines are the 3<r c.l. allowed experimental ranges, horizontal lines the resulting 
predictions for the fit (c2), see text. 



Acknowledgments 



This work was supported by Spanish grants FPA2005-01269 (MEC), ACOMP06/154 (Generalitat Va- 
lenciana), by Acciones Integradas HA-2007-0090 and by the European Commission Human Potential Pro- 
gram RTN network MRTN-CT-2004-503369. A.V. thanks the Generalitat Valenciana for financial support. 
W.P. is partially supported by the German Ministry of Education and Research (BMBF) under contract 
05HT6WWA and Fonds zur Forderung der wissenschaftlichen Forschung (FWF) of Austria, project. No. 
P18959-N16. 



Super-Kamiokande, Y. Fukuda et al, Phys. Rev. Lett. 81, 1562 (1998), hep-ex/9807003 



SNO, Q. R. Ahmad et al, Phys. Rev. Lett. 89, 011301 (2002), nucl-ex/0204008 
KamLAND, K. Eguchi et al, Phys. Rev . Lett. 90, 021802 (2003), hep-ex/0212021 . 
[MINOS Collaboration]. larXiv:0708l4 95 [hep-ex]. 
KamLAND Collaboration. larXiv:08"0T45 89 [hep-ex]. 

M. Maltoni, T. Schwet z, M. A. Tortola and J. W. F. Valle, New J. Phys. 6, 122 (2004); Online version 6 in 
arXiv:hep-ph/0405172 contains updated fits with data included up to Sep 2007. 
P. Minkowski, Phys. Lett. B67, 421 (1977). 

T. Yanagida, In Proceedings of the Workshop on the Baryon Number of the Universe and Unified Theories, 
Tsukuba, Japan, 13-14 Feb 1979. 

M. Gell-Mann, P. Ramond and R. Slansky, Print-80-0576 (CERN), Proceedings of the Workshop, Stony Brook, 
NY (North-Holland, Amsterdam, 1979). 

R. N. Mohapatra and G. Senjanovic, Phys. Rev. Lett. 44, 912 (1980). 
J. Schechter and J. W. F. Valle, Phys. Rev. D22, 2227 (1980). 

F. Deppisch, H. Pas , A. Redelbach, R. Ruckl and Y. Shimizu, Eur. Phys. J. C 28, 365 (2003) 
[arXiv:hep-ph/0206122j . 

G. A. Blair, W. Porod and P. M. Zerwas, Eur. Phys. J. C27, 263 (2003), lhep-ph /02 10058 . 
A. Freitas, W. Porod and P. M. Zerwas, Phys. Rev. D72, 115002 (2005), hep-ph/0509056 . 
M. R. Buckley and H. Murayama, Phys. Rev. Lett. 97, 231801 (2006) arXiv:hep-ph/0606088 . 
F. Deppisch, A. Freitas, W. Porod and P. M. Zerwas. larXiv:0712.036T1 jhep-phj . 

K. S. Babu, Phys. Lett. B203, 132 (1988). 



20 



D. Aristizabal Sierra, M. Hirsch and S. G. Kovalenko. larXiv:07 10.5699 [hep-ph]. 
L. J. Hall and M. Suzuki, Nucl. Phys. B231, 419 (1984). 

For a review on bilinear R-parity violation see, M. Hirsch and J. W. F. Valle, New J. Phys. 6, 76 (2004) 

[arXiv:hep-ph7 0405015 . 

W. Porod, M. Hirsch, J. Romao and J. W. F. Valle, Phys. Rev. D63, 115004 (200 1), [hep-ph/0011248] , 
D. Restrepo, W. Porod and J. W. F. Valle, Phys. Rev. D64, 055011 (2001), hep-ph/0104040 



M. Hirsch, W. Porod, J. C. Romao and J. W. F. Valle, Phys. Rev. D66, 095006 (2002), hep-ph/0207334 
A. Bartl, M. Hirsch, T. Kernreiter, W. Porod and J. W. F. Valle, JHEP 11, 005 (2003), lhep-ph/0306071 
M. Hirsch and W. Porod, Phys. Rev. D68, 115007 (2003), hep-ph/0307364 



D. N. Spergel et al [WMAP Collaboration], Astrophys. J. Suppl. 170, 377 (2007) arXiv:astro-ph/0603449 
S. Borgani, A. Masiero and M. Yamaguchi, Phys. Lett. B 386, 189 (1996) [arXiv:hep-ph/9605222 



F. Takayama and M. Yamaguchi, Phys. Lett. B485, 388 (2000), hep-ph/0005214 
M. Hirsch, W. Porod and D. Restrepo, JHEP 03, 062 (2005), hep-ph/0503059 . 
J. E. Kim, Phys. Rept. 150, 1 (1987). 

G. G. Raffelt, Stars as Laboratories for Fundamental Physics (University of Chicago, Chicago, 1996) p. 664. 
E. J. Chun and H. B. Kim, Phys. Rev. D60, 095006 (1999), [hep-ph/9906392] , 

M. Hirsch, PoS HEP2005, 343 (2006). 



M. Hirsch and W. Porod, Phys. Rev. D 74, 055003 (2006) [arXiv:h ep-ph/0606061 
C. S. Aulakh and R. N. Mohapatra, Phys. Lett. B119, 136 (1982). 
A. Masiero and J. W. F. Valle, Phys. Lett. B251, 273 (1990). 
Particle Data Group, W. M. Yao et al, J. Phys. G33, 1 (2006). 
R. Barbieri, S. Ferrara and C. A. Savoy, Phys. Lett. B119, 343 (1982). 

M. Hirsch, J. C . Romao, J. W. F. Valle and A. Villanova del Moral, Phys. Rev. D70, 073012 (2004), 
|hep-ph/0407269l . 

M. Hirsch, J~~C . Romao, J. W. F. Valle and A. Villanova del Moral, Phys. Rev. D73, 055007 (2006), 
|hep-ph/0512257] . 

F. Franks and S. Hesselbach, Phys. Lett. B 526, 370 (2002) [arXiv:hep-ph/0111285] . 

M. Hirsch, M. A . Diaz, W. Porod, J. C. Romao and J. W. F. Valle, Phys. Rev. D62, 113008 (2000), 
|hep-ph/0004"TT5] . 

J. C. RomaoTM . A. Diaz, M. Hirsch, W. Porod and J. W. F. Valle, Phys. Rev. D61, 071703 (2000), 
hep-ph/9907499| . 

M. A. Di'az,~KT Hirsch, W. Porod, J. C. Romao and J. W. F. Valle, Phys. Rev. D68, 013009 (2003), 



hep-ph/0302021 



'. F. Harrison, D. H. Perkins and W. G. Scott, Phys. Lett. B 530 (2002) 167 arXiv:hep-ph/0202074 . 
B. Mukhopadhyaya, S. Roy and F. Vissani, Phys. Lett. B 443, 191 (1998) [arXiv: hep-ph/980826 5^ 
S. Y. Choi, E. J. Chun, S. K. Kang and J. S. Lee, Phys. Rev. D 60, 075002 (1999) |arXiv:hep-ph/9903465|. 
F. de Campos, O. J. P. Eboli, M. B. Magro, W. Porod, D. Restrepo, M. Hirsch and J. W. F. Valle. larXiv:0712.2156l 

[hep-ph]. 

W. Porod, Comput. Phys. Commun. 153, 275 (2003) arXiv:hep-ph/030 lToT]. 
J. A. Aguilar-Saavedra et al, Eur. Phys. J. C 46, 43 (2006) |arXiv:hep-ph/0511344| . 

ATLAS collaboration, Rome meeting, 2004, https://twiki.cern.ch/twiki/bin/view/Atlas/RomeSUSYWiki 
B. C. Allanach, C. G. Lester, M. A. Parker and B. R. Webber, JHEP 0009 (2000) 004 [arXiv:hep-ph/0007009] . 
A. Barr, C. Lester and P. Stephens, J. P hys. G 29 (2003) 2343 [arXiv:hep-ph/0304226| . 
A. J. Barr, B. Gripaios and C. G. Lester. larXiv:0711.4 008 [hep-ph]. 
W. S. Cho, K. Choi, Y. G. Kim and C. B. Park. larXiv:071 1.4526 1 [hep-ph], 
A. J. Barr, private communications. 

As in the NMSSM, if the field $ is light and the coupling ho large, one has five neutralino states 
In the definition of the coefficient a given in [3!| there is a relative sign to the corresponding definition for the 
explicit case |42T ]. 
[59] We correct here a misprint in [34| . 



